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ERROR ANALYSIS OF BINARY RATE MULTIPLIER 
by George J. Moshos 
Lewis Research Center 


SUMMARY 

The binary rate multiplier is studied as a means of achieving approximate 
multiplication. The difference between the actual and desired output is defined 
as the error of the binary rate multiplier, and closed formulas are obtained 
for expressing this error in explicit form depending on the starting conditions 
of the binary rate multiplier counter. As a result of analyzing these error 
formulas, error bounds are obtained. 


INTRODUCTION 

An integral part of many special purpose digital computers used for real 
time control is the binary rate multiplier (ERM) (e.g., refs. 1 to 4). In 
these applications this unit is used as a means of scaling down a pulse stream 
to some specified fraction. A logic diagram of a BRM, which is built out of 
the standard logic elements shown in figure 1, is shown in figure 2(a). The 
NOR element shown in figure l(a) may have various number of inputs. The par- 
allel lines shown on one of the inputs of figure 1(c) are included to indicate 
that the AND circuit is intended to act as a pulse gate dependent on the level 
setting of the other line. The following brief description explains the opera- 
tions on the BRM. 

The input pulse stream is applied directly to the binary counter whose value 
is denoted by x n x n-l • • • X 2 X 1* Each flip-flop of the counter is operated as 

a trigger. For every two input pulses to a trigger, two output pulses are pro- 
duced; one pulse when the flip-flop makes a 0 to 1 transition, called an 

a pulse, and one when the flip-flop makes a 1 to 0 transition, called a (B pulse. 

The 3 pulse is used to trigger the next stage of the counter. The a pulses are 
gated through AND gates and mixed through a NOR element to produce the desired 
fraction of the input pulses. This simple mixing technique may be used because 
the a pulses from the various stages are separated in time from each other. 

This timing factor is shown in figure 2(b). 

This unit may be used to achieve approximate multiplication. In particu- 
lar, if Ax is the number of input pulses and y is a binary number less than 
1, Az, the number of output pulses, may be stated quantitatively as 


Az = y Ax 


(i) 



-8: 3D- 3- 


(a) Trigger flip-flop. (b) NOR. 

Figure 1. - Logic elements. 


(c) AND. 
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(b) Timing diagram. 

Figure 2. - Binary rate multiplier. 


The difference between the actual output of the BRM and Az, as given by equa- 
tion (l), is defined as the error of the unit. In this report, the formulation 
of error is explicitly given and studied. The analysis of the BRM is begun by 
deriving the approximate relation shown in equation (l). 


BINARY RATE MULTIPLIER 


The quantitative relation of a BRM may be expressed as follows : If Ax 

is the number of input pulses, the number of output pulses produced by the k^* 1 
stage of the counter is Ax • 2“^-. This multiplicative relation will remain 
valid over any interval for which Ax is a multiple of 2 k pulses. If y_k is 
the level setting of the k^ 1 stage AND gate, the number of output pulses that 


2 



may be gated through this stage will be y_& Ax • 2~ k . Since the output pulses 
from the various stages are simply mixed, the number of output pulses As of 
an n-stage BRM over any interval Ax, which is a multiple of 2 n pulses, will 
be the sum of all the pulses gated through all the stages. This output is 


n 




Az = Ax y) 
i=l 

The quantity y 
be written as 

n 

- E y-i*- 1 

1=1 

is a binary 



(1 

<< 

where the range 

of y is 

° < y < ! 


y_i2 _i ( 2 ) 

number. Therefore, equation (2) may 
Ax 

- 2~ n (3) 


in steps of 2" n . 


If y is constant over a Ax interval of 2 n pulses, the output given 

by equation (l) is exact. If y is constant over a Ax interval of less 

than 2 n pulses, however, this multiplicative relation may not be valid. In 

this latter case, the actual output depends not only on the values of y and 

Ax but also on the starting value of the BRM counter. If the output from a 

machine, whose BRM counter starting value is x s , is denoted by Az x , it can 

s 

be shown that the average output over all of the 2 n possible machines is also 
given by equation (l). This can be demo ns trated, as follows: If Az is the 

average output over all these machines, Az is, by definition. 


As = 




2 n 


( 4 ) 


2 n- l 

The sum ' As x given by equation (4) is the total pulse output over the 
3 

2 n different possible machines when each machine receives Ax input pulses 
and its AM) gates are set to the value y. It will be observed that the 
a transitions produced by these 2 n machines are the same as the ones produced 
by a single machine with 2 n Ax successive input pulses. For example, the a 
transition ending with counter value x is attained Ax times (once by each 
of the Ax machines whose starting value is prior to x in the counting se- 
quence) of the 2 n machines used in the average and also Ax times when 2 n Ax 
successive input pulses are applied to a single machine (since each counter 
value is traversed Ax times in this case). Therefore, the total pulse output 
over all the 2 n possible machines with each receiving Ax input pulses is 
equal to the pulse output of a single machine receiving 2 n Ax input pulses. 

In this latter case the pulse output is also given by equation (l) since the 
input pulse interval is a multiple of 2 n pulses. Therefore, 


3 



5 


( 5 ) 



0 2 4 6 8 

BRM counter value 

Figure 3. - Multiplication error resulting from y = . 101. 


on-i 



x s =0 


Combining equations (4) and (5) 
yields 

As = y Ax (6) 


MULTIPLICATION ERROR FORMULAS 

The error for the BRM may be 
defined as the actual output 
minus y Ax (i.e., the value pre- 
dicted by eq. (l)). Since the 
actual output changes only when 
the input pulses arrive, it is 
only necessary to consider the 
error at these discrete times. 

For example, figure 3 gives the 
output for a three-stage BRM to- 
gether with the values predicted 
by equation (l) when the BRM 
counter starting value is zero 
and y = .101. The difference 
between these two curves is the 


defined error. Yet, in discussing error, the error value will be given only 
when the abscissa is discrete values such as 0, 1, 2, etc. (which correspond to 
the BRM counter values). Moreover, it will be observed that when an output 
pulse is produced, the error may change by one quanta. Therefore, it is neces- 
sary to distinguish between the error immediately before the output pulse and 
the error immediately after the output pulse. When the BRM counter starting 
value is zero, these two values of the error will be denoted by F and E, 
respectively. When the BRM counter starting value is arbitrary, the starting 
value must also enter into the error formula as a parameter. The errors will 
be denoted prior to and after the arrival of the output pulse as H and G, 
respectively. It will be noted that the error defined this way makes E a 
special case of G, and F a special case of H. Nevertheless, this distinc- 
tion is maintained, since it is convenient for our subsequent discussion. 


Starting the BRM counter with zero gives the definition of the error E 
as the actual output after the output pulses are generated, minus y Ax. This 
difference, when only one stage of an n-stage BRM is gated by y, can be ex- 
pressed systematically in tabular form as shown in table I. 
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TABLE I. - MULTIPLICATION ERROR E OF BRM WHEN ONE STAGE IS GATED 


y-x 

a 1 

y-2 

= 1 

y -3 = 

1 

y-k = 

1 

x i 

E 

X 2 X Z 

E 

X 3 X 2 X 1 

E 

x k x k-l • • • X 3 X 2 X 1 

E 

0 

0 

0 0 

0 

0 0 0 

0 

00 ... 0 0 0 

0 

1 

1/2 

0 1 

-1/4 

0 0 1 

-1/8 

0 0 ... 0 0 1 

-1/2* 



1 0 

2/4 

0 10 

-2/8 

0 0 ... 0 1 0 

-2/2 k 



1 1 

1/4 

Oil 

-3/8 







10 0 

4/8 







10 1 

3/8 


• 





110 

2/8 

0 1 ... 1 1 1 

- (2 k - 1 - l)/2 k 





111 

1/8 

10 ... 0 0 0 

1/2 







10 ... 0 0 1 

(2 k -l - l)/2 k 







11 ... 1 1 1 

l/2 k 


An inspection of this table shows that the error associated with the various 
stages of a BRM may be expressed more concisely in algebraic form as shown in 
table II. 


TABLE II. - MULTIPLICATION ERROR E 
IN ALGEBRAIC FORM 


Stage 

E 

1 

y-i (x i/ 2) , , 

2 

y - 2 ^ 2/, 2 - x i/ 4 ) . 

3 

y-3( x 3/ 2 - x 2 / 4 - x n/ 8 ) 

k 

y_ k U k / 2 _ x k-i/ 4 ' • • • * x i/ 2k ) 


For an arbitrary value of y, the value of E is the linear combination 
of the values shown in table II. This bilinear form is shown in equation (7) 
for an n-stage BRM. The element subscripts of the Boolean vectors x and y 
(i.e., vectors whose elements are 0 or l), which are shown in this equation, 
correspond to the stage numbers of the BRM shown in figure 1: 



(7) 
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In the for mula tion of E, the maximum values of the output of the BRM were 
reflected at the points of discontinuities. It will be observed that just 
prior to these points the error is one quanta less than that shown by E. A 
formulation of F, in which the minimum values are reflected at the points of 
discontinuities, can be obtained in a manner similar to that for obtaining E. 
The quantity F, when only one stage of an n-stage BRM is gated by y, is 
shown in tabular form in table III (only a few cases are exhibited). 


TABI E III. - MULTIPLICATION ERROR F WHEN 
ONLY ONE STAGE IS GATED 


y_i = 1 y -2 = 1 

Cp Xl F CgC-^ x 2 x i 

0 0 0 0 0 0 0 
1 1 - 1/2 11 01 

10 10 
0 1 11 


y -3 = 


F 

c ; 

5 C ; 

> C 1 

X' 

5 X ; 

> X 1 

F 

0 

0 

0 

0 

0 

0 

0 

0 

- 1/4 

l 

1 

1 

0 

0 

1 

- 1/8 

- 2/4 

i 

1 

0 

0 

1 

0 

- 2/8 

1/4 

i 

0 

1 

0 

1 

1 

- 3/8 

l 

0 

0 

1 

0 

0 

- 4/8 


0 

1 

1 

1 

0 

1 

3/8 


0 

1 

0 

1 

1 

0 

2/8 


0 

0 

1 

1 

1 

1 

1/8 
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Comparing table III with table I shows that the values of F equal the 
values of E except at the points where the discontinuity occurs. At these 
points, F equals -l/ 2 , while the corresponding value of E equals +l/ 2 . The 
values of C shown in table III correspond to the 2 's complement of the x 
values. It will be observed that the values of F are identical in terms of 
C to the negative values of E. Consequently, it can be asserted that F in 
terms of C is just the negative of E: 


f « -(c^C;, . 



(9) 


An example will help clarify these formulas. In this example the value of 
y is . 101 , and the values of E and F are calculated for successive BRM 
counter values that are notated by the subscripts on E and F: 


/E n \ 

/000\ 

/ M 

/ 100 \ 

E 2 ' 


' 010 \ 

E 3 

= 

110 

E 4 


001 

E 5 J 


101 / 



V 011 / 
\ 111 / 

/’A 

/000\ 
/ 111 \ 

F 2 


f Oil \ 

F 3 

= - 

101 

E 4 


001 

E 5 j 


i 110 , 

E « / 

\ 010 / 

\ f ? / 

\ 100/ 


^1/2 -1/4 -1/8 
0 1/2 -1/4 

Vo 01/2 



/ 1/2 - 1/4 - 1/8 
I 0 1/2 -l /4 
V 0 01/2 




This example is shown in graphical form in figure 3 (p. 4 ). 


When a BRM counter starts out with an arbitrary value, the starting value 
must enter into the error formula as a parameter. These error formulas are 
given explicitly by equations ( 10 ) and (ll). In particular, these equations 
reflect the maximum (denoted by G) and minimum (denoted by H) values of the 
actual output at the points of discontinuities. In these formulations, x and 
xg represent the value and the initial value of the counter, and C and Cg 
represents the 2 *s complements of these values. The subscripts on those lite- 
rals represent, as before, the stage of the BRM. In equation (ll), xgR identi- 
fies the rightmost (i.e., lowest order) counter bit, whose value is 1; for 
example, for the counter value 100 , Xgj^y..^ = y_ 3, for counter value Oil, 

x sRy-R = y~i> etc: 
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1 



G = (x 1 - x S1 > x 2 ~ X S2 • • ♦ *n " x Sn )M 


( 10 ) 


7- 


n 


H - “ x SR y -R ~ ( C 1 " C S1' C 2 " C S2 * * 



( 11 ) 


The following two examples illustrate equations (10) and (ll) for y = .01 
and (xg^xgg ) = (0,l) and are calculated for successive BRM counter values that 
are notated by the subscripts on G and H: 




°o \( i/2 - 1 / 4 Y°) 

-1 J\o 1/2 Al/ 





'0 0\ /l/2 -l/4V0> 

1 -1 \ \0 1/2 Al ) 

1 0 / 



In these examples it will be noted that (xj_ - xgU may be 0, 1, or -1. 


MULTIPLICATION ERROR BOUNDS 

The maximum positive error and the minimum negative error for a n stage 
BRM whose counter starts out with zero may be obtained by an analysis of equa- 
tions (7) and (9), respectively. These values will then form a bound of the 
deviation of the BRM from that of exact multiplication. This analysis is pre- 
sented in appendix A. It is shown in that analysis that for an n-stage BRM 
these values are : 


E^ax^ 11 ) 
F min( n ) = 


= _7_ 

18 


+ S. + 


(-l) n 

n . on 


_ n _ (-l) n 

18 6 g . 2jn 


(12) 


(13) 


Equation (12) is plotted together with equation (13) in figure 4(a). As a 
by-product of developing equations (12) and (13), the points where these values 




Number of stages, n 


(a) Zero starting. (b) Arbitrary starting. 

Figure 4. - Multiplication error bounds. 


TABLE IV. - x AND y VALUES FOR ^ 


n 

Xy . . . X^ 

y_i . . . y_7 

X7 . . . x^_ 

y_i . . . y_7 

•^max 

2 

n 

11 

11 

11 

3/4 

3 

101 

101 

in 

m 

7/8 

4 

ion 

1101 

1101 

' ion 

17/16 

5 

10101 

10101 

non 

non 

39/32 

6 

loion 

110101 

noioi 

loion 

89/64 

7 

1010101 

1010101 

1101011 

noion 

199/128 


TABLE V. - x AND y VALUES FOR F . 


n 

Xy . . . XjL 

y_i • • • y_ 7 

X 7 • * • X-] 

1 

H 

1 

F min 

2 

01 

11 

01 

11 

-3/4 

3 

on 

101 

001 

m 

-7/8 

4 

0101 

noi 

0011 

1011 

-17/16 

5 

01011 

10101 

00101 

non 

-39/32 

6 

010101 

110101 

ooion 

loion 

-89/64 

7 

0101011 

1010101 

0010101 

noion 

-199/128 
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TABLE VI. - x, x s , AND y VALUES FOR 


n 

X S7 ’ ' ‘ ^1 

x 7 • • • X 1 

y-1 • • • y-7 

X S7 • • • *31 

x 7 • • • X 1 

y-1 • • • y-7 

^max 

2 

01 

10 

01 

00 

11 

11 

3/4 

3 

001 

110 

oil 

010 

101 

101 

9/8 

4 

0101 

1010 

0101 

0010 

1101 

ion 

23/16 

5 

00101 

11010 

01011 

01010 

10101 

10101 

57/32 

6 

010101 

101010 

010101 

001010 

110101 

101011 

135/64 

7 

0010101 

noioio 

0101011 

0101010 

1010101 

1010101 

313/128 




TABLE VII. - x, Xg y 

AHD y VALUES 

FOR H min 



n 

X S7 * ■ • X S1 

X 7 . . . x-L 

y-i • • • y-7 

X S7 • • • X S1 

Xy . . . Xjl 

y-i • • • y -7 

^min 

2 

n 

01 

n 

n 

01 

11 

-3/2 

3 

101 

on 

101 

in 

001 

in 

- 7/4 

4 

noi 

oon 

ion 

ion 

0101 

1101 

-17/8 

5 

10101 

01011 

10101 

non 

00101 

non 

-39/16 

6 

noioi 

001011 

loion 

loion 

010101 

noioi 

-89/32 

7 

1010101 

oioion 

1010101 

noion 

0010101 

noion 

-199/64 


occurred were found. These points are tabulated for various BRM in tables IV 
and V. 

Appendix B presents an analysis of a BRM whose counter starts out with an 
arbitrary value. The basis of this analysis is to use equation (10) to obtain 
the maximum positive error and to use equation (ll) to obtain the minimum nega- 
tive error. For an n-stage BRM these values are 

<W n ) = | + § - — (1 4 ) 

y 5 9 . 2 n 

®min( n ) = " g^^n-l ^ 15) 

These values form a bound for the generated round-off error. Fortunately, 
only two maximum and two minimum values may occur for a BRM (for n > 2), and 
therefore one can expect better results than would be predicted by these values . 
These values are plotted in figure 4(b) and are also presented together with the 
points at which they occur in tables VI and VII. 

The problem that has been considered in appendixes A and B and in this 
section is illustrated in figure 5. The actual and desired outputs for a two- 
and a three-stage BRM are plotted in these figures for all starting values . As 
is illustrated, finding the multiplication error bounds by graphical means is 
not trivial. The points labeled E max , F m i n , ( Wy , and H m -in in these simple 
cases agree with those predicted in appendixes A and B. 

The error formulas arrived at in appendixes A and B can be arrived at di- 
rectly by use of the exclusive OR operator. It was not used, however, because 
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y 1/8 p^ 4 rr^ i > 7 1 r° 


BRM counter, x 
(b) Concluded. 

Figure 5. - Concluded. 

some of the intermediate results obtained in the appendixes are interesting in 
their own right and would be bypassed by this alternate method of proof. 


CONCLUSIONS 

The BRM has been shown to achieve approximate multiplication. The error 
has been defined as the difference between the actual output and y Ax and is 
formulated as a bilinear expression. This error is shown to be dependent on the 
starting conditions of the BRM counter. The error formulas presented are ana- 
lyzed in detail and explicit error bounds are given. These bounds are shown to 
increase by approximately l/6 of a quanta per stage when the BRM counter start- 
ing value is zero, and by approximately l/3 of a quanta per stage when the BRM 
counter starting value is arbitrary. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 31, 1965 . 
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APPENDIX A 


MULTIPLICATION ERROR BOUNDS (ZERO STARTING) 

In this section the error equation of an n-stage BRM "whose counter starts 
with zero will be analyzed with the objective of obtaining tight error bounds. 
Nevertheless, some of the intermediate results that will be obtained in this 
section are interesting in their own right. Because this analysis is complex, 
formal methods will be used in this analysis. The basic outline is to use 
equation (7) to find the points where the maximum positive value is attained 
and then evaluate the equation at these points. In a similar manner, equa- 
tion (9) will be used to find the minimum negative value. 

The analysis is begun by stating and proving Lemma A.l. 

Lemma A.l : A sufficient condition for E given by equation (7) to attain its 

maximum value is that x^ = y_^. 

Equation (7) may be rewritten as a bilinear expression so that the terms which 
are dependent on either x^_ or y^ are grouped together. The quantity A in 
the resultant expression is independent of either xq or y_j_: 


E ( x i,y-i) = A + \ xj.y-i - \ Xi(f y-i- 


i-i + 7 y-i -2 



■ t y -i(t X i-1 + I Xi -2 + * * * + pri X l) (A3 

By direct evaluation, the value of this expression is as follows: 

E(0,0) = A 

K (i-°) = A - l(l y-1-1 + • • • + y-n) 

E(0,1) = A - |(| *i-l + • • • + -Zi x l) (At 

Moreover, for the specific case when i is n, the value of equation (Al) is 

E( 0,0) = A 
E(1,0) = A 


E(o,l) * A ' |(| ^n-1 + • ' • + ^ 

E(l,l) = A + f - + . . ■ 
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This Lemma is proved by observing that the value of E when Xj_ = y_i = 0 
is always greater than or equal to the value of E in both cases when Xj_ f y_ 
and that for the n^ 1 component the value of E is always greater when 
x n = y_ n = 1. Although it was not needed in the proof of this Lemma, it can be 
similarly shown that this is also true for the first component. 

Based on this Lemma, the maximum value of equation ( 7 ) will be obtained by 
finding the maximum of the quadratic form expression: 


Q(x-l,X2 . . . x n ) = (x-pXg . . . x n ) 

Theorem A.l : For all values of the components Xj_, 

Q(1,x 2 ,x 3 . . . Xn) > Q(0,x 2 ,x 3 . . . Xn) 
This theorem follows directly from the proof of Lemma A.l. 
Theorem A. 2 : For all values of the components Xj_, 

Q(x;l,x 2 ... Xn_i,l) > Q(xq_,x 2 ... Xn-i>0) 



This theorem follows directly from the proof of Lemma A.l. 
Theorem A. 3 : For all values of the components Xj_, 

Q(l,x 2 ,x 3 . . . Xn_}_,l) = Q(l,x 2 ,x 3 . . . Xn_i^l) 
where x^ is the complement of x^. 

The difference 

Q(l,x 2 ,x 3 . . . Xn_i^ l) ~ Q(l^ x 2- ,x 3 * * * x n-lU) = 

- (l,x 2 . . . x n _ 1 ,l)M 

may be written equivalently as 




(A4) 


(A5) 
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1 ^ 1 _ - 1 
2 " 2 2 " 2 n_1 


0 


1 _ JL L_ 

2 2 3 2 n ~ 2 


K = 



0 


1 


0 



1 _1_ 1 
2 " 2 1 " 2 n-i+1 


0 



Expanding equation (A6) and using the identity 


x i x j - x i x j = x i - x j 


to simplify the cross product terms give a typical term xj_ as 

( x i + - jr) 


(A7) 


(A8) 


Since (x-j_ + x-j.) = 1, equation (A6) is independent of the Xj_ variable, and the 
contribution from this term is 


1 JL. 
2 n-i+l - 2 i 


(A9) 


Similarly the contribution to the difference expressed in equation (A6) 
from the term involving x n _ i+ 2 _ is 


A _ 1 

2 i 2 n-i+l 


(A10) 


Consequently, the contribution to the difference expressed by equation (A6) by 
each element may be paired by the contribution from another element to cancel 
each other out of the expression. If n is odd, the middle term cannot be 
paired. But, since this is the (n + l)/z term, by equation (A9), its contribu- 
tion is 
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g(n+l)/2 " 2 (n+l)/2 

Therefore, the value of the difference shown by equation (A5) is equal to 
zero. This implies that 

Q(1,X 2 . . . x n _i,l) = Q(l,x 2 * * * x n _i,l) 

Lemma A , 2 : For v = xg^X 3 . . . Xj_ and a = 1,0, 1,0 . . • 1,0 

Q(l,v,0,a,l) > Q(l,v,l,a,l) 

where v and a are the component by component complement of v and a, 
respectively. 

By Theorem A. 3 it is noted that 

Q(l,v,l,a,l) = Q(l,v,0,a,l) 

Therefore, the difference between the two quadratic forms of the Lemma can be 
expressed as 


6 = Q(l,v,0,a,l) - Q(l,v,0,a,l) (All) 

Partitioning the M matrix of equation (All) so that Mp and M 2 are compat- 
ible with the vectors gives B as 

A 

5 = (l,v,0,a,l)Mpl v 

\0 


(l,v,0,a,l)M 2 A) 


- (l,v,0,a,l)Mp^A- (l,y,0,a,l)»%(i) (A12) 


But it will be noted that 


(l,v,0,a,l)M 



(l,v,0,a,l)M 



Therefore 


S = (l,7,0,a,l)M 2 Q _ (l ^^ 0 ,a,l)M 2 ^) (A13) 

It will now be proved by induction on the length of a that 5 > 0. It 
may be immediately verified that 5=0 for a of length zero. Assume that 
s k > where is the value of B when a is of length 2k. It will now 

be verified that &k+l > 
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1 

±2 

x 3 


x i 

°i+l 

1 i+2 

°i+3 


°i+2k-l 
x i+2k 
°i+2k+l 
1 i+2k+2 
°i+2k+3 
x i+2k+4/ 


-l/2 i+2k+ 4' 

-l/2i +2 k+3 

1 y 2 i+2k+2 


_-i /<?2k+5 

Jfftf k+4 

.l/ 2 2k+3 

_ 1 / 2 2k+2 









-l/2 i+2k+2 ' S 

. 1 / 2 i+2k+l 


.l/ 2 2k+3 

_-i /p2k+2 

-l/ 2 2k+l 

-l/2 2k 


(A14) 



The O's and l's in equation (A14) are subscripted to show their position 
in the vector, and the vector itself is displayed as a column rather than a row 
in order to show the correspondence between the terms that must be multiplied 
to form the value. 


Multiplying out the terms of equation (A14) yields 


3 k+l 


= 8b 


2 2k+3 ~ 2 i+2k+3 


_1 

2 2k+4 


i-2 

VA 

j=0 


x. 


1-0 


(A15) 


But, by direct evaluation 
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i-2 

1 1 - 1 12 1 
g2k+4 / Xi- ^ 2 2k+4 2 i_1 2 2k+4 2 2k+3 

j=0 


Therefore, &k+i > 0* 

Lemma A . 5 : For v = xg,xg ... xq and a = 1,0, 1,0 ... 1,0 

Q(l,v,l,0,a,l) > Q( 1, v, 0, 0,a, l) 


where v and a are defined as in Lemma A. 2. 


Proceeding in a manner similar to Lemma A. 2 shows by Theorem A. 3 that 

Q( 1, v, 0, 0, a, 1 ) = Q( 1, v, 1, 1, a, l) 


Therefore, the difference between the two quadratic forms of the Lemma can be 
expressed as 


S = Q( 1, v, 1, 0, a, 1 ) “ Q( 1, v, 1, 1, a, 1 ) (A16) 

Partitioning the M matrix of equation (A16) so that M]_ and Mg are 
compatible with the vectors gives 6 in the form 


& = (l,v,l,0,a,l)M^v^ + (l,v,l,0,a,l)Mg^ 


- (l,v,l,l,a,l)Mi^v^ - (l,v,l,l,a,l)Mg^a^ 


But expanding the first and third terms of equation (A17) shows that 

(l,v,l,0,a,l)M 1 ^ = (l,v,l,l,a,l)M 1 ^ 


Therefore, 


/ 0 \ _ / 1 > 

5 = (l,v,l,0,a,l)Mg[ a ] - ( l,v,l,l,a,l)Mg( a 

\l \V 


(A17) 


(A18 ) 


If & k is used to denote the value of 5 when a is of length 2k it will 
now be proved by induction on the length of a that 


8 k> 


- — x-i + —J=—r x i-l + • • ♦ + 


,2k+3 1 g2k+4 


ji+k+1 


+ 1 > 0 (A19) 


>i+2k+2 
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First, note that for a of length zero 



The notation used in equation (A20) is similar to that used in proving 
Lemma A. 2. Assume that 6^ > 0. It will now be shown that &k+l > 0: 



Multiplying equation (A2l) gives 
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5 k+l = & k “ 


(A22) 


1 - 


p 2k+5 X i “ 


1 - 

2 i+2k+3 X 2 


Using equation (A19) gives the right 

1 - 


> 


1 


22k+3 


2 i+2k+l x 2 + 2 i+2k+2 " 2 2k+5 


2 i+2k+4 

side of equation (A22) as 

1 - 


X-! - 


x 2 


> 


1 U 


Xi ”t . • . 


x 2 


> i+2k+4 


> 0 


2 i+2k+3 ° 2 i+2k+4 2 2k+5 - 1 - 2 i+2k+3 

Theorem A. 4 : There exists a v* such that for all v 

Q(l,v*,0,a,l) > Q(l,v,x i+1 ,a,l) 
where v and a are defined as before. 

First, note that Q(l,v**,0,a,l) > Q(l,v,l,a,l) because, suppose it were false, 
then there would exist a v** such that for all v 

* 


But from Lemma A. 2 


Q(l,v**,l,a,l) > Q(l,v,0,a,l) 


Q(l,v**,0,a,l) > Q(l,v**,l,a,l) 


Therefore, a contradiction exists. Moreover, note that 




Q(l,v**,0,a,l) > Q(l,v,0,a,l) 

that is, there is a largest. Therefore, Theorem A. 4 is proved by choosing 

either v** or v** for v*; that is, whichever makes Q(l,v*,0,a,l) the 
largest. 

Theorem A. 5 : There exists a v* such that for all v 

Q(l,v*,l,0,a,l) > Q(l,v,x i+1 ,0,a,l) 

First, note that Q(l,v**,l,0,a,l) > Q(l,v,0,0,a,l) because, suppose it were 

* 

false, then there would exist a v** such that for all v 


* 


Q(l,v**,0,0,a,l) > Q(l,v,l,0,a,l) 


But from Lemma A. 3 


Therefore , 


Therefore 


a contradiction exists. Moreover, a v® can he chosen so that 
Q(l,v**,0,0,a,l) > Q(l,v,0,0,a,l) 

Q(l,v*,l,0,a,l) > Q(l,v,x i+ i,0,a,l) 


Theorem A. 6 : There exists a v* such that for all v 

Q(l,v*,0,l) > Q(l,v,x 1+1 ,l) 

First, note that there exists a v** such that Q(l,v**,0,l) > Q(l,v,l,l), be- 

It* 

cause; suppose it were false, then there would exist a v such that for all v 

Q(l,v**,l,l) > Q(l,v,0,l) 

But by Theorem A. 3 

Q(l,v**,0,l) = Q(l,y**,l,l) 

** 

Therefore, a contradiction exists. Moreover, a v** can be chosen such that 

** 

Q(l,v**,0,l) > Q(l,v,0,l) 


Therefore, 


Q(l,v*,0,l) > Q(l,v,x i+1 ,l) 

It will now be demonstrated by an example how these theorems can be used 
to obtain the value of x so that the error is the maximum positive value. 
Consider a seven-stage BUM. By Theorem A.l and Theorem A. 2 

Q(l,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,l) > Q(x 1 ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,x 7 ) 


by Theorem A. 6 

Q( 1 ,x|,X 3 ,X 4 ,X 5 , 0 , 1 ) > Q( l,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,l) 

by Theorem A. 5 


Q(l,x£*,xf*,x|*, 1,0,1 > Q(l,x|,x|,x|,x§,0,l) 


by Theorem A. 4 




by Theorem A. 5 


** 


* 


Q(l,x£*,l,0,l,0,l) > Q(l,xf%x§*,0,l,0,l) 
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by Theorem A. 6 


Q(l, 0,1,0, 1,0,1) > Q(l,x**, 1,0, 1,0,1) 

Putting these inequalities together results in 

Q(l,0,l,0,l,0,l) > Q(x 1 ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,x 7 ) 

Moreover, by Theorem A. 3 

Q( 1, 1, 0, 1, 0, 1, 1 ) > Q( X x , x 2 , x 3 , x 4 , x 5 , Xg, x ? ) 

Using the theorems in the pattern illustrated by the example makes it 
easy to verify that the maximum positive value will occur at the points shown 
in table IV ( p . 9 ) . 


The maximum positive value of the error may be expressed concisely as 
follows: Let E may (k) denote this value for a k-stage BRM. If k is odd. 




max v 


x 3 




Ik 

°k+i . 
x k+2/ 


-i/a 

vh j 


(A23) 


Evaluating equation (A23) yields the difference equation 

E max (k+2) = E max (k) + — ^1 + ^+2^ 

Solving equation (A24) for an n-stage BRM in terms of E max ( 1 ) yields 


(A24) 


Ei 


'max 


(n) - E max (l) - i + S. - — — 


->n 


But E max (l) — l/2. 

Therefore, the maximum positive error of an n-stage BRM, where 


E max( n ) 



n is odd, is : 
(A25) 


If k is even. 
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1 1 

3 gk+2 


( A26 ) 


Solving this difference equation for an n-stage BRM in terms of E max (2) and 
then evaluating the resultant expression for ^max^) = 3/4 yield 


•^max( n ) 


J~ + £ + — 1 ■ 

18 6 9 . 2 n 


(A27) 


Combining equations (A25) and (A27) gives a elosed-form equation for the maxi- 
mum positive error of an n-stage BRM: 




E, 


'max 


(n) ’ + | + A3f 

18 6 9 • 2 n 


( A28 ) 


The minimum negative value for an n-stage BRM can be obtained by applying 
equation (9). Comparing the form of equation (7) with equation (9) shows that 
the previous results can be utilized with a slight modification. In particular, 
the value of the minimum is equal to the negative of the maximum and occurs at 
points that are the 2 j s complement of the maximum value. Consequently, the 
minimum negative values will occur at the points shown in table V (p. 9). 
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APPENDIX B 


MULTIPLICATION ERROR BOUNDS (ARBITRARY STARTING) 

The error formulas given by equations (10) and (ll) express the multiplica- 
tion error of a BRM whose counter starts with an arbitrary value. Equation 
(10) is the error formula resulting when the maximum value of the actual output 
is considered at the points of discontinuities. Equation (ll) is the companion 
equation resulting when the minimum value of the actual output is considered at 
these points. In this section, these error formulas will be analyzed with the 
objective of obtaining error bounds for a BRM with this added degree of freedom. 
First, equation (10) is analyzed to obtain the maximum positive error of an 
n-stage BRM. 


It is convenient for this discussion to define a vector b so that 



= M 



(Bl) 


n 

Theorem B.l : For all x^ and Xg k in equation (10), G < |b_jj. Moreover, 

i«l 

n _ 

G = ^ jb_.j, where x^ denotes the complement of Xc 


if y-k = = x sk> 


i=l 


Sk* 


The elements of the vector defined by equation (Bl) are 


b -i - 1 y-i - (3 y-s + 1 y-3 + • • • - ^ y-n) 

b -2 - \ y-2 - (I y-3 + 1 y-4 + • • • + y-n) 

b -k - 1 y-k - (f y-k-i + 1 y-k-2 + • • • - 53+1 y -») 

B_ n = 2 y~n 

Since l/2 > l/4 + . . . + l/2 n . 
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< 0 

if y_k = 0 


> 0 

if y_ k = 1 

( B2 ) 


It is next noted that each element (x^ - xg^) may have only three possible 
values; that is , 0,1, or -1. This may be verified by direct computation: 


x k 

x Sk 

x k - x Sk 

0 

0 

0 

0 

1 

-1 

1 

0 

1 

1 

1 

0 


Since |xfc - x S k| < 1 then 


G - (x x - x sl )b_ 1 + (x 2 - x S2 )b_ 2 + 


+ 


~ x Sr>-n < 


n 


E l b _i| 


(B3) 


(B4) 


A sufficient condition for G to attain the upper bound of equation (334), 
n 

that is, ^2 |b_jj, is that 
i=l 


( x k - x Sk) = -1 

if b_ k < 0 

= +1 

if b_ k > 0 


Combining equations (B2), (B3), and (B5) results in 


y-k 

b -k 

x k - x Sk 

x k 

x Sk 

0 

£0 

-1 

0 

1 

1 

>0 

+1 

1 

0 


__ n 

Therefore, y_^ = x^ = Xg^. is a sufficient condition for G = ^ |b_jj • 

i=l 


As a consequence of Theorem B.l the maximum of G, denoted by G y , is 
n 

such that 0 mx = max ^ |^-i| # Procedure to be followed is to find the 

y i=i 


n 

value y where the maximum of |^-i| a *ttained and then evaluating this 

i=l 

function. In order to aid this analysis the notation b_^(y) is introduced, 
where y = * * * y_ n an< ^ b_^(y) denotes the value b_^ for the vector 


25 




II || || I 1 1 llllll IHIIIIIHI Hill I ■■llllll II II II nil III mil IIHMI ■■Mill I •■■Ml 


(y_l,y_2 • • • y_ n )* For an n-stage BRM, all possible b_i(y) values may be 
obtained by multiplying M with all possible values of y. This particular 
matrix is called B n . 



Vl(D 

b o(l) 


b_i( 2 ) 
b .(2) 


V-n (1) 


b_ n (2) 


b .-^ 11 - 1 )' 


'A 


b. 2 (2 n - 1) 


(B6) 


b-nt 2 ” 



A few examples will help clarify equation (B6) 
Two- stage BRM: 



Three-stage BRM: 



V° 0 V 


( r 


0 0 0 1 


0 110 


yi 0 1 


1 1 l\ 
Oil 
10 1/ 



1 _3 

4 8 

1 1 

2 4 


1 3 

2 8 


1 

4 


1 

2 


0 I 



Four -stage BRM: 



_4 _5_ __6_ __7_ 1 _7_ _6_ _5_ _4_ _3_ _2_ Jl\ 

16 ~16 "l6 16 2 16 16 16 16 16 16 16 \ 


1 

2 


3 

8 


1 

4 


1 

8 


0 


B 


4 “ 



Be 


0 

0 


1 

4 

1 

2 


1 

2 

0 


1 

4 

1 

2 


0 

0 


B 3 
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Theorem B,2 : For all values of y, |k_i(y)| = |b_.j_(2 n -y)|. 

This result follows by induction on the number of stages k, It was shown as 

true for the k = n - 1 case 


1 2 n-1 - 1 1 \ 
rr * ♦ • \ 

2 2 n 2 n 

0 

• B n-1 

6 / 

where this case is partitioned to show its structure . This theorem is obviously 
true for the first row. The b_j_(y) element for the n - 1 case is now the 
b_i_i(y) and the b_j__i(2 n + y) elements of the n case; and the b_j_(2 n - y) 
element of the n - 1 case is now the b_i_]_(2 n - y) and the b_i_i( 2 n + 2 n - y ) 
elements of the B n case. By the induction hypothesis |b_]_(y)| = |b_i(2 n -y) 
for the n - 1 case. Therefore, these elements for the n bb case yield 

|b_i_i(y)| = |b_ i _ 1 (2 n+1 - y) | (B7) 

and 

|*>_i_i(2 n - y)| = |b-i-i(2 n + y) | (B8) 

Substituting u = 2 n - y into equation (338) gives 

|b_i_i( u )| = |'b_ i _ 1 (2 n+1 - u)| 
which completes the proof. 

Lemma B.l : There exists a y* in the domain 

010, . , 00<y*<0111. . *11 so that 


£| b -i(y*)| > D | b -i(y)| 

i i 

This theorem states that G-joax is attained in the domain 

010, , ,00<y<011, , ,1L As a result of Theorem B.2, the search 
for a point where G max is attained can be immediately restricted to the y 
domain 0 < y < 1 0 0 . . , 00, Consider B n for these values of y: 


an example for the k = 2 case. Assume it i 
Then the k = n case is 


/_ 1 _ 2 - 2 n_1 - 1 
/ — — » • • 

2 n 2 n 2 n 


B n - 


V 


B. 


n-1 
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I s ®n-l* Because of Theorem B.2, the absolute values of the elements in 
are identical to the absolute values of the elements E|r. Moreover, since the 
elements of first row, that is, |b_]_(y)|, increase as y increases, then for 



0 

each column sum to the left of . there is a column sum to the right that 


\ 0 / 

exceeds it. Therefore, G^y must lie to the right. 

As a result of Theorem B.2 and Lemma B.l, there must be at least two values 
of y where G mPV is attained. For an n-stage BRM, the y value corre- 
sponding to G m Py on the left of y = 1000 ... 00 will be called the 

one on the right of y = 100 ... 00 will be called R n . 

Lemma B.2 : For an n-stage BRM 

0R n _ 1 <L n <011. . .11 

This result follows from the proof of Lemma B.l. Since the first row; that is. 
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t>_i(y)|^ is increasing, G mY must lie between the right maximum of B n _^ and 
;he rightmost value of Lemma B.l. 


Lemma B.5: For an n- stage BRM 


0100. . . 0 < Ln < 0 1 Ln_ 2 

Consider Bj^ for the values of y of Lemma B.l; that is, 
0100. . • 0 < y < 0 1 1 . . .1: 



g n-2 + 2 
„n 



B 


n 


1 2 n ~ 2 - 1 2 n ~ 2 - 2 
2 2 n “ 2 n ”^" 



B n-2 


0 



The column by column sums of the absolute value of the elements of the 
first two rows are 


3 2 n “l + 2 n ” 2 - 1 2 n “l + 2 n " 2 - 2 

2 n ' ' 2 n ' ‘ ‘ * 

Therefore, this is a decreasing sequence. Since by Theorem B.2 and Lemma B.l, 
B n _2 attains a maximum for at least two values, then G mav must lie between 
the leftmost value of Lemma B.l and 01L n _2. 

Theorem B.5 : For an n- stage BRM 

0R n-l < L n < 01L n-2 


This theorem is the combination of Lemma B.2 and Lemma B.3. 

Theorem B.4 : R n and are unique and 0R n ^ 1 = L n = 011^2 • This theorem 

follows immediately from the proofs of Lemma B.2 and Lemma B.3 by using the 
principle of strong induction as the method of proof; that is, assume it is 
true for k < n and prove for the case k = n. The values of y where 
is attained can be obtained by using Theorem B.4. These values are listed in 
table VIII. 

The BRM counter value and starting value corresponding to the y values 
listed in table VIII can be obtained by Theorem B.l. These values are listed 
in table VI ( p . 10 ) . 

An equation for G my as a function of n may be obtained by a procedure 
similar to that used to obtain E mx . In particular, if the pattern established 
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IS ATTAINED 


TABLE VIII. - VALDES OF y WHERE G maY 


OH'ra -2 

0 B n-l 

n 

y-i • • • y-6 

y-i * • ■ 

1 

2 

3 

4 

5 

6 

on 

0101 

oion 

010101 

0X1 

0101 

01011 

010101 


h 


y-i ■ * • y-6 

y-i ■ * • y-6 

i 

i 

01 

n 

on 

101 

0101 

1011 

01011 

10101 

010101 

101011 


in table VIII is used a difference equation may be written for n even, and a 
difference equation may be written for n odd. Combining the solutions of 
these difference equations gives G^y as a function of n: 


max 



(-D n 

9 • 2 n 


(B9) 


Noting the similarity between the rightmost term of equation (ll) to that of 
equation (10), establishes immediately a minimum error bound when the BRM 
counter starts with an arbitrary value: 


B min( n ) 1 ~ ~ 


n , (-l) n 

3 9 • 2 n 


(BIO) 


A tight error bound may be obtained as follows: Expanding equation (ll) 

gives the equation for H in the form 

H = ~Xg R y_ R - £(Cq - Ggj_)b_q + (Cg - Cgg)b _2 + . . . + (C n - Cg n )b_ n J (Bll) 


In the initial value of the BRM counter, xg R identifies the rightmost 1. 
It may be simply argued that, if a binary number has a rightmost 1 in position 
R, its 2’s complement also has a 1 in that position and, moreover, has zeros 

at all positions j, where j < R. Therefore, xg R = Cg R and Cgj = 0 for all 

j < R. Equation (Bll) may be expressed as 

-H = C gR y_ R + ^C 1 b_ 1 + C 2 b_ 2 + . . . + 0 R _jb_ R+ ^ + (C R - Cg R )b_ R 

+ [( C R+1 “ c SR+l) b -R-l + • • • + ( c n " c Sn) b -n] ( B12 ) 

or alternately as 

-H = C-jb.q + C 2 b_ 2 + . . . + C R b_ R + (C R+1 - Cg R +i)b.. R _ 1 

+ . . . + (C n - Cg n )b_ n + C SR (| y_ R + i y_ R -l +•••+) (B13) 
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Therefore 


< | c l b -l| + | c 2 b -2| + • • • + | C R b -R | + |( C R+1 " °SR+l) b -R-l| + • • • 

+ | (C n “ C Sn^ b -nl + C Sr(| ^-R + \ ^-R-l + • • •) ( B14 ) 

The upper bound of -H is attained when equality is attained in equa- 
tion ( B14 ) • It can be demonstrated by an argument similar to that used in 
Theorem B.l that the conditions for equality are 


y -j c j 


for all j 


and 


Cgj = Cj for all j > R 
Therefore , equation (B14) may be written as 


l 


C2 


-H = (C 1} C Z ,C 5 


C n )M 


\ °n / 


- ( 0,0,0 . . . 0,Cp+2_,Cj}+2 


C n )M 



C Sr(2 °R + I °R+1 + ’ • *) (B15) 


The last two terms in equation (B15) are always nonnegative. Moreover, the 
sum of these two terms is nondecreasing as R decreases. Since by Theorem A.l, 
C]_ = 1 is a condition for maximizing the first term of equation (B15), it must 
also be a condition for maximizing equation (B15) itself. 

Equation (Bll) may be rewritten with this condition as follows: 

“ H - y_i + ( c 2 - C S2> C 3 “ C S3 * 4 * C n " C Sn) M l • 1 (B16) 



Therefore, the equation for H m i n can be immediately written as: 


HminU) = -1 - CW(n-l) = -ig - - -LDl 


9 . 2 


,n-l 


(B17) 


Based on this analysis Cg^, xg]_, Cq_, xj_, and y_]_ are equal to 1 for 
Hmin* ^e remaining stages are determined to maximize G for an (n - l)-stage 
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TABLE IX. - C, y, AMD C s WHERE 
H pin IS ATTAINED 


n 

y-i - • * y-6 

Cg . . . c 1 

C S6 • * ’ °S1 

2 

xi 

n 

01 

3 

101 

101 

on 

4 

ion 

noi 

oon 

5 

10101 

10101 

01011 

6 

loion 

noioi 

ooion 

n 

y-i • * * y-6 

c 6 • • • C 1 

C S6 ' ' ' C S1 

2 

n 

11 

01 

3 

in 

m 

001 

4 

noi 

ion 

0101 

5 

non 

non 

00101 

6 

110101 

ioio n 

010101 


BRM. The values of y, C, and Cg, where H m j_ n is attained, are listed in 
table IX. The BRM counter value and the initial value for these %in values 
are the 2's complement of the preceding numbers. These values are tabulated 
in table VII. 
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